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Abstract 

We discuss tlie Overhauser effect (particle-hole pairing) versus the BCS effect (particle- 
particle or hole-hole pairing) in QCD at large quark density. In weak coupling and to 
leading logarithm accuracy, the pairing energies can be estimated exactly. For a small 



q_) , number of colors, the BCS effect overtakes the Overhauser effect, while for a large num- 

i-^ ' ber of colors the opposite takes place, in agreement with a recent renormalization group 

argument. In strong coupling with large pairing energies, the Overhauser effect may be 

^^ . dominant for any number of colors, suggesting that QCD may crystallize into an insu- 

3 , lator at a few times nuclear matter density, a situation reminiscent of dense Skyrmions. 

The Overhauser effect is dominant in QCD in 1-1-1 dimensions, although susceptible to 

quantum effects. It is sensitive to temperature in all dimensions. 
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1. Introduction 



Quantum chromodynamics (QCD) at high density, relevant to the physics of the early 
universe, compact stars and relativistic heavy ion collisions, is presently attracting a renewed 
attention from both nuclear and particle theorists. Following an early suggestion by Bailin 
and Love [0, it was recently stressed that at large quark density, diquarks could condense 
into a color superconductor |||], with potentially interesting and novel phenomena such as 
color-flavor locking, chiral symmetry breaking, parity violation, color-flavor anomalies, and 
superqualitons. 

At large density, quarks at the edge of the Fermi surface interact weakly thanks to 
asymptotic freedom. However, the high degeneracy of the Fermi surface causes perturbation 
theory to fail. As a result, particles can pair and condense at the edge of the Fermi surface 
leading to energy gaps. Particle-particle and hole-hole pairing (BCS effect) have been 
extensively studied recently [||, ^ . Particle-hole pairing at the opposite edges of the Fermi 
surface (Overhauser effect) Q has received little attention with the exception of an early 
variational study by Deryagin, Grigoriev and Rubakov for a large number of colors Q, and 
a recent renormalization group argument in O. The scattering amplitude between a pair 
of particles at the opposite edges of the Fermi surface peaks in the forward direction, a 
situation reminiscent of the forward enhancement in Compton and Bhabha scattering. 

In retrospect, it is surprising that the Overhauser effect in QCD has attracted so little 
attention. In fact, the Schwinger model [Q shows that when a uniform external charge den- 
sity is applied, the electrons respond by screening the external charge and inducing a charge 
density wave, a situation analogous to a Wigner crystal 0, ^, ^. Similar considerations 
apply to QCD in 1-1-1 dimensions H. In 3-1-1 dimensions, dense Skyrmion calculations with 
realistic chiral parameters yield a 3-dimensional Wigner-type crystal with half-Skyrmion 



symmetry at few times nuclear matter density [10, 11 1. At these densities, Fermi motion is 
expected to be overtaken by the classical interaction |jl2| . A close inspection of these results 
shows the occurrence of scalar-isoscalar, pseudoscalar-isovector and vector-isoscalar charge 
density waves in an ensemble of dense Skyrmions. 

In this paper we will show that in dense QCD, the equations that drive the particle- 
hole instability at the opposite edge of a Fermi surface resemble those that drive the particle- 
particle or hole- hole instability in the scalar-isoscalar channel, modulo phase-space factors. 
In section 2 we motivate and derive a Wilsonian action around the Fermi surface. In section 
3 we obtain expressions for the energy densities and pertinent gaps in the 0"*" channel with 
screening, thereby generalizing the original results in Q. In section 4, we analyze the 
decoupled equations for large chemical potential without screening. The effects of screening 
for arbitrary Nc as well as temperature are discussed in section 5, in overall agreement with 
a recent renormalization group argument Q. In section 6, we discuss the Overhauser effect 
in QCD in lower dimensions. Our conclusions and suggestions are given in section 7. 



2. Effective Action at the Fermi Surface 



To compare the Overhauser effect to the BCS effect, we wih construct a Wilsonian 
effective action by integrating out the quark modes around the Fermi surface, in the presence 



of smooth bifocal fields. An alternative would be the quantum action [13|. At large chemical 
potential, most of the Fermi surface is Pauli-blocked, so the quasiparticle content of the 
theory is well described by such an action. Incidentally, our analysis should provide a 
useful alternative to a brute-force lattice QCD analysis. Indeed, an effective formulation of 
lattice QCD along the lines of the heavy-quark formalism is possible and will be discussed 
elsewhere [p^]. 

The starting point in our analysis is the appropriate QCD action in Euclidean space 
with massless quarks 



d^x 



\iF^uf + i^il.d^ - 74/i)V' - iJ^K 



(1) 



and the colored current 



^^ = 5V'7myV'. (2) 

In Euclidean space, our conventions are such that the 7-matrices are hermitean with 
{l^l■,lv} = '^^ij.v For sufficiently large ^, we will assume g^Nc <^ 1. We have omitted 
gauge-fixing terms and ghost-fields. In what follows, we will analyze (|l]) in the one-loop 
approximation with the gluon field in the Feynman gauge. The approximation, as we shall 
show below, is equivalent to the resummation of the ladder graphs in the particle-particle 
or particle-hole graphs. The effects of screening will be dealt with by minimally modifying 
the gluon propagator, ignoring for simplicity vertex corrections as in [Q]. The issue of gauge 
fixing dependence will be briefly discussed at the end. 

In the one- loop approximation with screened gluons, the induced action is 

^V = Y / d^^ dS J^{x)V^,ix - y)J^iy) + j d^x^; d^^^^l^ , (3) 

where 5^ = {di, 82,83,84 — 11). The screened gluon propagator {V^^} = {Ve, T>m) is 

^.,M(.-.) = /i0^,^e-(-). (4) 

Perturbative arguments give m\/{gii)'^ = m'j^/{gfj,)^ ps Nf/27r'^ and m'j^j/m'j^ « 7r|g4|/|4q|, 
where mo is the Debye mass, niM is the magnetic screening generated by Landau damping 
and Nf the number of flavors |15| ^. Nonperturbative arguments suggest m\,m\.[ — > 



mz 



^jq^ 1 16] where for simplicity, the difference between electric and magnetic channels is 



ignored. We expect Kqcd ^ "n^* < itt-e in the case Nc = 3, as lattice simulations for the 
*^ Throughout we will refer to ttim abusively as the magnetic screening mass. 



gluon propagator at finite /x are not yet available. We note that the perturbative screening 
vanishes at large Nc- 

To proceed further with ^ we need to Fierz rearrange the J J term in (^). This is 
equivalent to summing ladder graphs with relevant quantum numbers. Specifically, 

o 

+ g^J2^0'[Hx)Mo'riy)]'0{x-y)[riy)Mo'^ix)] (5) 
O' 

with Co = —1/9 and Cc = +1/36 for the operators 

[i^{x)M.cij'{y)\ = i^a,a,i{^){l5)af3eiheijCi^lh^j{y), (6) 

respectively, with Nj = N^ = 3. These quantities involve matrices active in color (a, 5, • • •), 
flavor (z,j, •••) and Dirac space (a, /?, •••). Mo is the vertex generator for particle-hole 
pairing in the 0"*" channel (i.e., Overhauser), while Mc is the vertex generator for particle- 
particle and hole-hole pairing in the color-flavor locked (CFL) channel (i.e., BCS). Only 
these two operators will be retained below, unless specified otherwise. The gluon-propagator 
in matter is 

V{x-y) = ^VE{x-y) + ^VM{x-y) . (7) 

The weightings follow from minimal substitution in matter with 2 electric and 2 magnetic 
modes. We note that the present Fierzing is particular, since it selects solely the Ic in 
the qq channel and the 3c in the qq channel |17|. For arbitrary Nc > 3 and Nj > 2, the 



coefficients — | and ^ become, respectively, —^(1 — ■^) ■ -^ and 2^ • ^ • minfjv N ) ' where 
the single factors refer, in turn, to the results of the color Fierzing, the flavor Fierzing and, 
of course only for the second expression, the Fierzing related to color-flavor locking ^. To 



compare to the more conventional decompositions through 3c x 3c = Ic + 8c for qq and 

in - M.J- and ^^^+1 . 1 . 1 



3c X 3c = 3c + 6c for qq, with respective weights — A(l — 1^) ■ 7^ and — s±l • i • 3_ 



we introduce also the vertex generator [g] 

[i^ix) Mb V''(y)] = ^Pa,aAx) (75)a/3 (A2)a6(T2)»jCV^ftbj(2/) ■ (8) 

We note that (0) does not lock color and flavor as it stands; a color-flavor locking as 



described in footnote #2 still has to be performed, such that flnally the corresponding 



*^At least partially, even if Nf ^ Nc, a locking can be achieved by Fierzing the antisymmetric tensor 
in color times the corresponding one in flavor into the tensor 'M.ai,bj = SaiSbj — SajSti with the pertinent 
weight l/inm{Nc,Nf) in the combined color-flavor space. The latter operator has ^n{n— 1) eigenvalues -1-1, 
^n{n -I- 1) — 1 eigenvalues —1, one eigenvalue n — 1, and Nc x Nf — v? eigenvalues 0, where n = rma{Nc,Nf). 
Thus, in the BCS case, the fermion determinant (h6|) acquires the color-flavor weight 2n(n — 1), while the 
corresponding value in the Overhauser case is the standard NcNf factor. 



coefficient becomes Cb = {Nc + I) / {8NcTain{Nf , N^)} *^. This brings about the important 
issue of whether Fierzing is a unique operation on 4-fermi interactions. The answer is 
no [^, ^. This nonuniqueness would of course not be important if an ah-order calculation 
were to be performed for any Fierzing set, but is of course relevant for truncated calculations 
as is the case in general. Each Fierzing corresponds to summing a specific class of ladder 
diagrams in the energy density, see e.g. [^, |23| ]. 

Introducing a hermitian bilocal field T,{x,y) and a non-hermitian bilocal field T(x,y), 
we may linearize the Fierzed form of the J J term by using the Hubbard-Stratonovich 
transformation, e.g. 



with 



exp ( fg / d xd y ['0(x)V'(y)] V{x - y) [il){y)il){x)] 

(iS(x, y) exp f -5's - / d'^x (fy i;{x) T,{x, y) ip{y) j (9) 

S^ = ^Jdxdy,^^-—^ (10) 

and similarly for F. As a result, the action in the quark fields is linear and the functional 
integration can be performed. The result is the following effective action for the bilocal 
fields 

S = Ss + Sr- iTr InF , (11) 

where 

p ^ / {7 • 5 - M74} S{x -y) + MoS(x, y) iV^x, y)C^Mc \ 

[ iC^T{x,y)Mc {7-a + ^74}<5(x-2/) + MoS(x,y) j' ^ 



The factor of 1/2 in (|ll|) is due to the occurrence of ijj and ijj^ through the Fierzing into 1^ 



and 3c |17]. This renders naturally the Gorkov formalism applicable to the present problem 
even at /i = 0. Note that Mo = Ic x If x 1_D and Mc = Sq x Sp xj^, with the subscripts 
C, F, D short for color, flavor and Dirac. We should stress that the effective action ( |ll|) 
is general. The third term is the Hartree contribution of the quarks to the ground state 
energy at large chemical potential, while the first two terms remove the double counting in 
the potential (i.e., Fock terms). 

To analyze the Overhauser and BCS effects in parallel, we make simplifying ansatze 
for the bilocal auxiliary fields. Since the unscreened gluon interaction in both cases peak in 
the forward direction, we may choose 



T,{x,y) = 2 cos 



p , Xfj. + Un 



a{x — y) = 2 cos 



P.(^)]/i|^e-.-(-.)F(,) 



*^In fact, expression (H) is the operator considered in Ref. |lq, vM where color and flavor are uncoupled 
and only the two-flavor case is considered (see also pi EOl). Thus the corresponding coeflicient is just 
{Nc + l)/4Nc, since the flavor- Fierzing factor can be ignored, as it eventually cancels against a corresponding 
factor resulting from the fermion determinant. Note, furthermore, that our color-flavor coupling scheme is 
difi'erent from the one recently introduced in Ref. |2l| for arbitrary numbers of flavors. 



T{x,y) 



2 cos 



P„ 



Vf^ 



g{x — y) = 2 cos 



P„ 



Vfi 



d^q 



(2vr)^ 



iq-{x-y) 



G{q) 



(13) 



where P^ = (Pf,0) and |P_p| = 2^. Y'p points in the original direction of one of the 
quark. F{q) and G{q) are even functions, F[q) is real, since S(x,y) = S(y, x)*, and G{q) is 
complex, since T^{x,y) = T{y,x)* . The relative momentum q satisfies \q\ < \P/2\ = fi. The 
bilocal field F characterizes a BCS pair of zero total momentum. S characterizes a wave of 
total momentum 2/x. This is the optimal choice for the momentum of the standing wave for 
which the holes contribute coherently to the wave formation. As a result the gap opens up 
at the Fermi surface, with n as the divide between particles and holes. In both cases, the 
pairing involves a particle and/or hole at the opposite sides of the Fermi surface. Indeed, 
in terms of (O) the linear terms in the bilocal fields are 



d xd yip{x)T,{x,y)il^{y) 



Va 



(see Ref.i]) and 



d\ 
(2vr)4 



V' 



^+q) F{q)i^{^+q)+^(X+q) F{q)^(-^+q 



(14) 



\ d'^xd'^y V''(a;)«75r(x,?/)V'(y) +^(a:)F'''(x,y)i75^''(?/) 



^V, 



{2nr 



l.T 



P 



P 



fP 



+ V 



T 



P 



r [ -^-qj Cij5Giq)i^i^-+qj +iji^- + q^ iG* {q)-i^C^I^' (^---g 
P 



P 



Cil^Giq) ^ 



+q]+i^ 



^+q) iG*{qh5G^^(^-q 



(15) 



m 



that 



where V4 is the 4-volume. 

Following [H], we introduce fermion fields ip{-^P/2 + q) and ^'^(ibP/2 — 
are independent integration variables in the relevant region of the momentum Ig^l <C |P|/2. 
Hence, the quark contribution around the Fermi surface can be integrated. The result is |25| 



/ 



det F = exp 



Vi Tr In 



V 



-iQ+,^a^ Fj^q) iG*{q)mc 

F{q) -iQ-,^a^ -iG*{q)Mc 

-iG{q)Mc -^Q;,^^m F{-q) 

iGiq) Mc F{-q) -iQ*.,f,a^ 



\ 



, (16) 



where Q± = ±^ + q and Q± = (Q± , Q^. —ifJ^)- For each entry in momentum space 
q, the determinant in ( p!6[) is over an {8 ■ Nc ■ Nj) x [8 ■ N^ ■ A^j)-matrix. The matrices 
cr^ = (i(T, 1) and a^ = {—ia, 1) are defined in terms of the usual Pauli matrices a. The 
detailed analysis of the coupled problem (16) with the full Fermion determinant will be 
discussed elsewhere |25|. 



**Note that we define tp'^ik) = C'ip'^{k) in terms of the Euclidean charge conjugation operator C — 7472, 
whereas in Ref. |24] ip^ik) = CijF {—k). As usual, T stands for "transposed". 
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3. Gap Equations 



A qualitative understanding of the Overhauser effect versus the BCS effect can be 
achieved by studying the phases separately, and then comparing their energy densities at 
large quark density. Setting G = yields, for the Overhauser pairing, an energy density 



9^4 g^J V{x) 



(2^ 



In 



QlQl + 2F2 Q+Q_ + F^ 



S, 



pot,E 



+ 5] 



kin,S 



(17) 



which is in agreement with the result derived originally in Q. Setting F = yields, for the 
BCS pairing, an energy density 



Sr 1 [,.^\9{xr 
36^4 g^ J V{x) 



dS 

(2vr)4 



In 



QlQ 



*2 



^pot.r + ^kin.r 
(18) 



which is similar to (p^). In writing the last equation we have assumed that \G{q it P) 
|G(ib/i)| K, 0. The gap equations for both cases follow by variation. The result is 



F{p) = 2g^ 



d^q 



(2vr)4 

for the Overhauser gap, and 

d^- 



V{p - q) 



2F{q)[Q+Q_+F^{q) 
QlQl+2F^{q)Q+Q^+F\q) 



(19) 



G{p) 



si/ 



i2^Y 



V{p - q) 



2G(g)(Q+Q; + |G(g)|- 

QlQf + 2\G{q)\'^ Q+Q\ + \G{q)\^ 



(20) 



for the BCS gap. If we were to use the antisymmetric vertex operator (|8|) then we would 
have 25^/3 — > 4^(^/3. For the latter, we have checked that the results (p!9|-pO|) agree with the 
Bethe-Salpeter derivation in the ladder approximation to order p,^ . In our notations, the 
leading order effects are of order p., the next to leading order effects are of order ^j? and the 
next-to-next to leading order effects are of order //~^. For the screened gluon propagator 
we have the alternatives 

1 1 1 



V{q) 



"^ q^ + m\ 



+ k 



"^q^+mlj 



1 



"^ q^ + im 



+ i — 



1 



imt 



(21) 



for the perturbative and nonperturbative assignments respectively. 

The present construction is valid for an arbitrary number of colors with or without 

screening, thereby generalizing the original analysis in Q. The outcome can be analyzed 

variationally, numerically or even analytically to leading logarithm accuracy. Using the 

following momentum decomposition around the fixed Fermi momentum P at the Fermi 

surface, 

P q P 

«ii = ^r' ^^=q-^ii^' (22) 
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and assuming that the relevant values of the amplitudes of the bilocal fields are small (i.e., 
F,\G\ <C ^), we may further simplify the kinetic part in the energy densities eqs.(p!7|-p^). 
Specifically, 



^kin.S 






In 



2 ~|2 



,^ + F2(,) + {,4 + t} 



The simplified gap equations are 
and 




^(^)-i^V(S^^(^-^) 



Ql + ilW + iY 



G{q) 

ql + \G{qW + {q\\ + f^y 



(23) 



(24) 



(25) 



For both pairings, the simplified energy densities 5e r at their respective extrema are 



5e 



\G{q)\' 



2 f d^q /I 



;|G|a|G|-l In 1 + 



36^4 3J (27r)4 V2'-'-'-' J ' y ' ^2 + (^^ + |l,2 

We now proceed to evaluate F, G to leading logarithm accuracy. 



4. Unscreened Case: Large Nc 



In this section we consider the gap equations (^^]25|) in the absence of screening. 
In the perturbative regime, we note that mE,M ~ l/-^c) and this approximation may be 
somehow justified in large Nc B. Hence, 



F{p) « 2g^ I 



d\ 1 

(2^)4 {p-qf 



F{q) 



ql + F^{q) + {q^ + ^~^^ 



2i 



and 



2 2/- dU 1 

^(^) " 3^ i (2^ (^^3^ 



Gjq) 

ql + \G{q)\' + {q\\ + ir 



(26) 



(27) 



For the Overhauser pairing, if we assume the propagator to be static, the q^ integra- 
tion can be performed by a contour-integration with the constraint that 



q±\'<2f,€g = 2fiJqf.+F^{q^^). (28) 



Hence 



,2 /- ... ^(l\\'l,.J, , 2fM, 



Fip,)«,'J^ *llif >"(l+(^^j P») 



2 

with h'^ = j-^. In general, we have 



^ - 87r2 Iv N, 



8^2 V Af2 



for Fierzing with Mc and M^ ^ respectively. Eq. (|29| ) is essentially a one-dimensional 
'fish-diagram' with logarithmically running couplings. This feature is preserved by screening 
as we will show below, in agreement with the recent renormalization group analysis in Q. 
Following [^^, the resulting equations are readily solved by defining the logarithmic scales 
X = ln(2/i/p||), y = ln(2/x/g||), xq = ln(2/x/Fo), and rewriting 

F{x) « h? (2x r dy F{y) - H dy yF{y) + f dy yF{y)] . (31) 

V Jx Jx Jo / 

Since F"{x) = -2h^F{x) with F(xo) = --F(O), then F{x) = -Focos{V2hx) g, |ol. The 
coefficient Fq follows from F'{xq) = 0, with ^/2hxQ = vr. Hence Fq = F{xo) and xq = 
ln(2/x/Fo) = Tr/{y/2h). Thus 

f„.2.exp{--^} (32) 

which is exactly the result established in Q using the Mc Fierzing and elaborate variational 
arguments. Note that the pairing energy Fq <C A_|_ <C /U follows from an exponentially small 
region in transverse momentum ( [2^ ) as required by momentum conservation, see Fig. lb. 
Typically A_l = \/2fj.Fo as originally suggested in Q. 

For the BCS pairing, the transverse momentum is not restricted as shown in Fig. la. 
This is best illustrated by noting that the BCS equation in (|27| ) can be further simplified 
through the following substitution 



q2 



^^11 + ^ 



^+^ 



// . (33) 



This amounts to taking into account the effects of curvature around the fixed Fermi mo- 
mentum P/2 defined by the standing wave. The trade (33) allows for a larger covering of 



the Fermi surface, although for A_l = 2/x the terms that are dropped are only subleading 
for qn <^ q'j_. We have checked that this substitution does not not affect our analysis in 



the leading logarithm approximation. Shifting momenta to Q = q + P/2 and K = p + P/2 
yields 



2„2 f d'^Q ^,r. ^^ r G{Q-P/2) 



GiK-P/2)^-,^jAQ^iK-Q) 



Ql + \G{Q-P/2)\^ + {\Q\-f,y 



(34) 



For a constant gap, the Q-integration diverges logarithmically. As most of the physics 



follows from |Q| = /i, this divergence can be regulated |18|], with no effect on the leading- 
logarithm estimate of the pairing energy. Hence, 



with eg = ^/q?, + |G(g||)p following from the contour integration over Q4. The prefactor 
reads h^ = (p' JVIti^ ^ and in general ^ 

2 _ 5W 2 \ min(Arf,Ar,)-l 



87r2 \NJ 2 

corresponding to Fierzing with Mc and M^ respectively. Notice the similarity between 
( p^ and (|35[), especially in the one-dimensional reduction of the equations. In terms of the 
logarithmic scales, the BCS equation reads |2C] 

G{x) « 2hl (^x j'J dy G{y) + j^ dy yG{y)^ . (37) 

Since G" {x) = -2hlG{x) with G(0) = 0, then G{x) = Gosm{V2Kx). The coefficient 
Go follows from G'{xq) = with ^/2h^XQ = 7r/2. Hence Gq = G{xq) and, because of 
xq = ln(2/i/Go), 

Note that Go is enhanced relative to Fq, if Nc = 3. They both become comparable for 
A^c > 4 in the Mc-Fierzing case with the Overhauser effect dominating at large Nc ^, as 
originally suggested in 0. 



We note that the i in ( P4[) (Overhauser) versus no i in (|25| ) (BCS) stems from the 
kinematical difference between the two pairings, hence a difference in the phase-space in- 
tegration due to momentum conservation as shown in Fig. la and lb. In weak coupling, 
both gaps are exponentially small. The energy budget can be assessed by noting that the 
phase space volumes are of order: /i^ Go (BCS) and ^Fq (Overhauser). Hence, the energy 
densities are 

|«-/.f3, f-VG^ (39) 



*^In |20| and footnote 1 of |5| color and flavor are uncoupled. Hence h?, — ^-7(1 + 7^)- In fact, this value 
can also be reproduced by the Fierzing with Ms for the special case Nf = 3. 
*®In the Ms-Fierzing case, the Overhauser effect only dominates in the large Nc limit, if Nf < \Nc. 
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(a) 



(b) 




X 



Figure 1: (a) fraction of the Fermi surface used in BCS pairing; (b) fraction of the Fermi 
surface used in the Overhauser pairing with one standing wave; (c) fractions of the Fermi 
surface used in the Overhauser pairing with two orthogonal standing waves. 



In weak couphng, the BCS phase is energetically favored up to N^ ~ 10 in the unscreened 
case and for one standing wave for the Mc Fierzing. Under M^ Fierzing, we have an 
additional constraint on the number of flavors, e.g., A'^^ < ^Nc for large Nc- Remember 
that further nestings of the Fermi surface by particle-hole pairing are still possible as shown 
in Fig. Ic, causing a further reduction in S's/V4. A total nesting of the Fermi surface will 
bring about 47r/i^/A^ « ji/Fq patches, hence Sy^/V^ « —^j^Fq. The BCS phase becomes 
comparable to the Overhauser phase for N^. ^ ^ (see, however, footnote #6). Finally, we 
note that in strong coupling, both gaps are a fraction of /i. 



5. Screened Case: Finite A',- 



In the presence of electric and magnetic screening, which are important in matter, the 
situation changes significantly. While the original variational arguments in Q were tailored 
for the unscreened case, our formulation which reproduces exactly their unscreened results in 
the leading logarithm approximation, generalizes naturally to the screened perturbative and 
nonperturbative cases in a minimal way. Indeed, using (|2^Pq) and the pertinent transverse 
cutoffs, we obtain for perturbative screening, 



h f°° 



X In 



^(^ii: 



lqf\+FHqn) 



1 + 



A2 



{p\\-q\\y + m' 
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X In M 1 + 



G{q\\) 
i + \G{q\\)? 

A2 



{p\\-%f 



1 + 



A3 



m. 



\P\\-<1\\? + lmjj\P\\-Q\\\ 



(40) 



and for nonperturbative screening 







G{p\\) 



dq\\ 



G{q\\) 



In 



1 + 



A 



In 



1 + 



-q\\f + iml 
A2 



{P\\ -^||)^ + 



imi 



(41) 



where the transverse cutoffs are A_|_ = y^2^ue^ (Overhauser) and A_i_ = 2/i (BCS) respectively. 
The cutoffs are exactly fixed in weak coupling, and reflect on the fractions of the Fermi 
surface used in the pairing. 

In the BCS case, the transverse cutoff is large. Hence A_|_ > ruE^'mM and the loga- 
rithm in (|40|) may not be expanded. Dropping 1, we obtain to leading logarithm accuracy. 



Gn 



4A'^ \ _VM 



vrmj 



(42) 



The results for the BCS gap is the same as the one reached in ||2^, ^, 18 1 ^ if we were 
to Fierz with M^ instead of Mc*. Note that (|4^ ) is smaller than ( |38|) as expected. For 
nonperturbative screening, the result is 



Go ~ A±e 



r{ln{l+A\/rrii)} ^ 



(43) 



with A_L/m* = 2/x/m=K » 1. 

In the Overhauser case, the transverse cutoff is reduced in comparison to the BCS 
case due to momentum conservation for fixed 3-momentum for the standing wave. The 
equation can be rearranged into the form 



/l2 






^^(911) 



In 



2/ieg 



M\-q\\y 



5h^ 



+ -yo''^II^^K(^^^ 



rrii 



(44) 



where we have approximated irm'j^/A ~ m'^ and used static, but perturbatively screened 
propagators. The effects of Landau damping through the magnetic gluons result into an 

*^Modulo the dimensionless constant feg in llSl. 
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unscreened interaction but with a reduced strength h? -^ /i^/6. Eq. (^) can be solved to 
leading logarithm accuracy using the logarithmic scales as defined above. Specifically, for 

X < Xm = ln(2^/m£;), we get 

/ pX r PX r pX \ 

F{x)^h?[2xj dyF{y)-J dyyF{y) + J^ dyyF{y)\ (45) 

as in the unscreened case with xi = xq, and for x > Xm 

F{x) « — i^2x j^ dy F{y) - j^ dy yF{y) + j^ dy yF{y)j + C . (46) 

Here xj^jj = ln(2/i/i<j^jj) and the constant C is given by 

c=^r,,„£Mi„f 'n- ,y (47) 

6 Jo eg \max[qf,,m%)l 



The solution to (|45|-[4q) is 



(x) = 


= FLCOs{^/2h{x - Xl)) 


for 


^ ^ ^m •> 


(x) = 


= FRCOs{h{x - xr)/V3) 


for 


X > ^m • 



(48) 

We note that for x < Xm or p^^ > niE, screening can be ignored to leading logarithm 
accuracy and xl = vr/v^/i as before. For x > Xm or pn < niE, screening cannot be ignored 
to leading logarithm accuracy. Continuity at Xm fixes xr, so that 

p^ ^ p^ cos{a/2/.(x^-x^)} ^ ^_,/^, ^o^^^^^fe)-^} (49) 

'' cos {h (x„ - xr)/V3} cos {h [in (^) - In (|^)] /^s} ' 



Eq. (49) defines a transcendental equation for Fr/2^ as a function of Nf,Nc and h (equiv- 



alently /x), i.e. 






where we have used mE/^/j, ~ h J Nj /Nc , with 

— r PS Rs — In . ^ H In In . ^ (51) 

to two loops. For fixed Nc and in weak coupling (/i ^ 0), there is no solution to ( |50[ ) as 
can be seen by inspection. This corresponds to a screening mass with power suppression, 
e.g. rriEJ \i. ~ h. However, a solution can be found in weak coupling but large A'^c, when 
approximatly 

for which Fr ~ Fe- Through N^^ this corresponds to a screening mass with exponential 
suppression, e.g. mEJ [i ~ ^-■^N'^h^ 

13 



To assess the minimal value of Nc for which there is a solution to ( |50D , it is useful to 
note that the solution ( p8| ) is invariant under the shift 

with similar shifts in the scales Xm,L,R, implying the existence of a family of solutions that 
depend parametrically on Xm,L,R and Aj^. The harmonic equation satisfied by F{x) is scale 
invariant, hence of the renormalization group type ^, The scale xr is fixed in terms of 



XL,m by demanding that the logarithmic derivatives of (|48| ) (with pertinent shifts) match 
at Xm- Thus 

-j=i&n\^—r={xm-XR)\=iMi[y2h{xrn-XLyj , (54) 

with 



Xij = — — j arctan VQ tan (V2h{xm — xl)) +mod7r>+Xm- (55) 

The lower bound on Nc or equivalently the upper bound on the electric mass follows from 

THE = A^e-^™ = [ ^ J (^) e-^™ < 2/i (^^ e^^'™ = 2/ie-^'ir^'" , (56) 

where A_l and An are now exponentially small scales characterizing the spread in p± and 
P||. The inequality in (|5^) follows from the geometrical constraint An > A^/2;U discussed 



above (see (Pq) and also Fig. lb). Up to the rescaling (53), the maximum An for which 
there is a solution ( ^8|) with positive semi-definite gap, corresponds to F{x\\) = 0, i.e. 
a^ll = xr + \/3TT/2h. (The alternative solution X|| = xl + \/27r/4/i does not generate a 
maximum bound.) After inserting the latter and (55) into c = v2^niin(x|| + Xm), we 
determine the minimum as c ~ 2.5051 and the lower bound for Nc (upper bound for the 
electric mass niE) as 

Nf 



This result is in overall agreement with a recent renormalization group estimate [0] ^. In 
particular, for /x = SAqcd, we find Nc > SSANf. 



The case of nonperturbative screening can be addressed similarly by noting that ( [44[ ) 
is now 

2 Jo eg \{p\\ -q\\r + mtj 



** Indeed, f{x) = -F'{x)/F{x) satisfies f (x) = f\x) + 2h'^ for x < Xm and f'{x) = f{x) + h^ /i for 
X > Xm, which are the renormahzation group equations derived in [|g], after the identification h -^ h/^/2. 
A similar observation extends to the BCS case, where g{x) — —G'{x)/G{x) satisfies g'{x) = g'^{x) + 2hl 
(unscreened) and g'{x) — g^{x) + h'i/3 (screened), in agreement with the renormahzation group equations 
derived in [ po[ . 

*® After the identification h -^ h/\/2 and /j, ^ 2fi in the prefactor of tue in H. 
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For X < Xm or p|| > m^, the screening in (|5q ) is inactive. Hence F{x) = —FLCos{\/2hx), 
while for x > Xm or p\\ < m^, the screening overwhelms the leading logarithm accuracy with 
F{x) = const. Continuity at Xm requires that Xm = xl. Hence 

m* = 2/u e ^h (59) 

which is the maximum tolerated nonperturbative screening mass for an Overhauser pairing 
to take place. 

Finally, we can qualitatively analyze the effects of temperature on the Overhauser 
effect by considering the distribution of quasiparticles at the Fermi surface. At finite tem- 
perature r, the pairing energy becomes 



-<..)=» t/:^..^-(o^)-(* 



{p\\-<i\\f 



+ ^ L ^^11 ^ '- i (pii-,i)^;^.(T)^ j *-^ {^) ' (^°) 



with the temperature dependent screening mass IH 



m%{T) = m%+[N, + ^]^-—. (61) 



Even at large Nc the screening mass is finite. We conclude that at finite temperature, the 
Overhauser pairing is rapidly depleted by screening for any value of N^. 



6. Pairing in Lower Dimensions 

The results we have derived depend on the number of dimensions. Indeed, the QCD 
analysis we have carried out when applied to 1+1 dimensions yield the following energy 
gaps 

F{q) 1 

/O 



F{p) ^ h^ dq 

Jo 



oo 



{p - qY + ml ' 

with the replacement (/^/Stt^ — > g^ /2t: in /i^ and h^. Remember that F{q) and G{q) have 
been defined as even functions. In deriving ( |62[ ) we have followed the same logic as in 3+1 
dimensions, thereby ignoring self-energy insertion on the quark line, and the gauge-fixing 
dependence on the gluon propagator. While these two effects cancel in color singlet states 
(Overhauser) p6|, they usually do not in color-non-singlet states (BCS) except for the case 
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of A'^c = 2 |22]. In 1+1 dimensions g^ JIt^ has mass dimension, and there is only electric 
screening with m|; ~ -^/<7^ In (/^/s)- Clearly, 

Fo f« A e""!/'^' > Go ~ A e"'"!/'^* . (63) 

The dominance of the Overhauser effect over the BCS effect whatever N^^ stems from the 
fact that the Fermi sm'face reduces to 2 points (±/x) in 1+1 dimensions, with no phase space 
reduction for the former. Since both the Overhauser and BCS phase break spontaneously 
chiral symmetry at finite density, the existence of the Overhauser phase may rely ultimatly 
on large N^- The Overhauser effect is dominant in the Schwinger model where Gq = 
because of the repulsive character of the Coulomb interaction ^^, confirming the results 



in [0, §, ^. The case of QCD in 2+1 dimensions will be discussed elsewhere. 



7. Conclusions 



We have constructed a Wilsonian effective action for various scalar-isoscalar excita- 
tions around the Fermi surface. Our analysis in the decoupled mode shows that in weak- 
coupling, the Overhauser effect can overtake the BCS effect only at large N^ in the scalar- 
isoscalar channel, in agreement with a recent renormalization group result [^. The BCS 
pairing is more robust to screening than the Overhauser pairing in weak coupling. The BCS 
analysis was carried out for both the CFL and the antisymmetric arrangements for arbitrary 
^c ^ 3, Nj > 2, ignoring the superconducting penetration lengths since the electric and 
magnetic screening lengths are smaller than the London and Pippard lengths (for type-I 
superconductors). In strong coupling, the Overhauser effect appears to be comparable to 
the BCS effect, especially if multiple standing waves are used, allowing for further cooper- 
ative pairing between adjacent patches. This is particularly relevant for pairings with large 
energy gaps which are expected to take place at a few times nuclear matter density p!9[| . 

Our effective action is better suited to the use of variational approximations as dis- 
cussed in y], and leads naturally to exact integral equations by variations, especially in 
the presence of interactions with retardation and screening. It would be interesting to re- 
peat our analysis at nonasymptotic densities using instanton-generated vertices to address 
the Overhauser effect. Indeed, for instantons the cutoff is fixed from the onset by their 
inverse size. As we have shown here, the Overhauser pairing, much like the BCS pairing 
by magnetic forces [pO|], relies on scattering between pairs in the forward direction that is 
kinematically suppressed in the transverse directions (in fact exponentially suppressed Q). 
Since the instanton interaction is nearly uniform over the Fermi sphere, we expect a geomet- 
rical enhancement in the BCS pairing in comparison to the Overhauser pairing. We recall 



#10 



In the Schwinger model m% = g'^ /2n independently of /i. 
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that in the latter the interaction is enhanced by a factor of order Nc- Which one dominates 
at a few times nuclear matter density and A'^c = 3 is not clear a priori. Instantons in the 
vacuum crystallize for Nc > 20 ||2^ in the quenched approximation, and 3 < A'^c < 20 in 
the unquenched case. The crystallization is likely to be favored by finite fj, as the quarks 
are forced to line-up along the forward X4-direction. 

It is amusing to note that the crystal phase breaks color, flavor, and translational sym- 
metry spontaneously, with the occurrence of color and flavor density waves. In many ways, 
this situation resembles the one encountered with dense Skyrmions [|lO| (strong coupling), 
suggesting the possibility of a smooth transition. In the process, color and flavor, respec- 
tively, may get misaligned [p8|, resulting into color- flavor- locked charge density waves in a 
normal (large gaps) phase. The Skyrmion crystal at low density may smoothly transmute to 
a qualiton crystal at intermediate densities, with crystalline structure commensurate with 
the number of patches on the Fermi surface. We note that the crystalline structure in 3 -|- 1 
dimensions may only show up as rapid variations in the response functions at momentum 
2fi. This is not the case in 1-1-1 and 2-1-1 dimensions [^. 

Although we have carried out the analysis using Feynman gauge with minimal changes 
for the electric and magnetic screening, we expect our estimates of the gap energies to be re- 
liable since a close inspection of the equations we derived when reinterpreted in Minkowski 
space, shows that the quoted results originate from the forward scattering amplitude of 
quarks around the Fermi surface. The latter is infrared sensitive in the unscreened case and 



gauge independent, the exception being in 1-1-1 dimension |26, 22|. The similarity between 
forward particle-particle and particle-hole scattering resembles the similarity between for- 
ward Compton and Bhabha scattering. This is what makes 2^ and opposite sides to the 
Fermi surface so special between a particle and a hole. 

Finally, it is amusing to note that following either the Overhauser or BCS pairing, 
the quark eigenvalues of the QCD Dirac operator would suggest a novel rearrangement that 
is characterized by novel spectral sum rules. They will be reported elsewhere. Our use 
of the effective action at the Fermi surface is more than a convenience for the study of 
QCD at large quark chemical potential. Indeed, given the shortcomings faced by impor- 
tant samplings in lattice Monte Carlo simulations at finite quark chemical potential, and 
also given the importance of Pauli blocking for the non-surface modes, we believe that a 
convenient formulation of QCD on the lattice should make use of Fermionic fields projected 
onto the Fermi surface, much like the ones used in the present work, and in the spirit of the 
heavy-quark formalism []T4|. 
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